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It has been found recently that propagators, e.g. the cross-correlation spectra of the cosmic 
fields with the initial density field, decay exponentially at large-fc in an Eulerian description of the 
dynamics. We explore here similar quantities defined for a Lagrangian space description. We find 
that propagators in Lagrangian space do not exhibit the same properties: they are found not to 
be monotonic functions of time, and to track back the linear growth rate at late time (but with a 
renormalized amplitude) . These results have been obtained with a novel method which we describe 
alongside. It allows the formal resummation of the same set of diagrams as those that led to the 
known results in Eulerian space. We provide a tentative explanation for the marked differences seen 
between the Eulerian and the Lagrangian cases, and we point out the role played by the vorticity 
degrees of freedom that are specific to the Lagrangian formalism. This provides us with new insights 
into the late-time behavior of the propagators. 

PACS numbers: 98.80.-k, 98.80.Bp, 98.65.-r 



I. INTRODUCTION 

Although the details of the evolution of the large-scale 
structure of the universe are probably affected by the 
presence of baryonic matter, in the context of a dark 
matter dominated cosmic fluid, it is thought that the 
global properties of the matter distribution at cosmolog- 
ical scales are essentially determined by those of a self- 
gravitating dust fluid. 

A complete understanding of the development of grav- 
itational instabilities in such a fluid is still however an 
open problem. It is one of the central issues for the 
study of structure formation in observational cosmology 
and this is for instance what pure TV-body cosmologi- 
cal simulations attempt to solve. The Vlasov equation, 
that is the fluid limit of the Boltzmann equation, en- 
tirely describes this system (see [1] or [2] for details). 
This equation of motion applies to the so-called Eulerian 
description of the dynamics, where the fluid properties 
are described through functions of fixed space-time co- 
ordinates (such as density and velocity fields). However, 
there exists an alternative description where the system 
is defined by the trajectories of particles, labelled as a 
function of their initial positions. This is the Lagrangian 
formalism, which takes advantage of the particle descrip- 
tion of the fluid. While this may not be a convenient 
description of the dynamics in the fully developed non- 
linear regime, it gives good insights of the dynamics in 
the early stages of the development of the gravitational 
clustering. The widely used Zel'dovich approximation [3[ 
corresponds for instance to a description of the displace- 
ment field based on its linear approximation. 



There exists a standard perturbative approach to study 
the development of gravitational instabilities beyond the 
linear approximation. This approach, and the main re- 
sults it led to, is described to a large extent in [2|. While 
it can can be useful for some specific observables, it fails 
to provide effective tools for describing the evolution of 
quantities such as the density power spectrum beyond the 
linear regime. Then, one still needs to use semi-analytic 
prescriptions. The ones that are mostly used now, e.g. 
the so-called Peacock and Dodds formula [H or the Smith 
et al. formula [H , originate either from the near universal 
transform advocated in ft|. or are based on an even more 
empirical construction, the halo model (see Q)- It is to 
be noted though that these prescriptions offer predictions 
for the power spectrum with relatively low accuracy, at 
the level of 10%, and are insecure in cases of non stan- 
dard cosmological models. Clearly there is a need to do 
better! 

Recently there has been a revival of Perturbation The- 
ory techniques (see [1, H, E3, EH and also [l2j for an 
overview of these ideas). In particular the Renormalized 
Perturbation Theory (hereafter RPT) formalism intro- 
duced in [§[ suggests a new scheme for the construction 
of perturbation theory expansions. It has been success- 
fully applied to the shape of the two-point propagator, 
and consequently to the two-point density power 
spectrum One of the core objects of this approach 
are the so-called propagators. They can be viewed as the 
cross-correlation between the cosmic fluids (that can ei- 
ther be the local density contrast or the peculiar velocity 
divergence) and the initial density field. In particular, it 
has been found that these correlators decay exponentially 
in the large-fc limit (where k is the Fourier mode of in- 
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terest). This result has been obtained analytically from 
a partial resummation of diagrams - in a perturbation 
theory point of view - that are thought to be the lead- 
ing contributors of the high-fc behavior of this quantity. 
It has been furthermore confirmed in numerical simula- 
tions. 

The aim of this paper is to consider similar quantities 
in the Lagrangian description of the dynamics. Thus, our 
goal is not to reconstruct the real spa ce p ower spectrum 
from Lagrangian variables (as done in [la ]) but to extend, 
to other objects of interest that arise in the Lagrangian 
framework, exact PT results. We first recall in section HT1 
the basic ingredients of this description. To compute the 
high-fc limit of the propagators we then assume Gaussian 
initial conditions and that the same set of diagrams will 
provide us with the leading contributions. To be more 
specific, those diagrams are those in which all loop s are 
connected to the principal line. As shown in [161 ]. and 
explained in details here, this approximation amounts to 
linearize the motion equation for a mode evolution while 
the low-fc modes act as a random stochastic background. 
As we show in section Hill for the 2D dynamics and in sec- 
tion [IV] for the 3D, although the modes of this stochastic 
background - assumed to be of Gaussian statistics - are 
in infinite number, their effects can be recast as those 
of a finite number of Gaussian random variables. This 
method reveals extremely powerful. We explicitly show 
the results it leads to for the 2D and 3D-Lagrangian prop- 
agators. We summarize in the last section what we have 
learned from these calculations. 



II. LAGRANGIAN APPROACH 

A. Equations of motion 

In Lagrangian approaches the global properties of the 
fluid are reconstructed from the individual particle tra- 
jectories, x(q, t), labeled by their initial Lagrangian co- 
ordinate q. Thus, the Eulerian comoving position x at 
time t reads as 



*(q, t), 



(i) 



where *(q, t) is the displacement field. Note that in 
Eq.(JTJ we use the property that in standard cosmolog- 
ical scenarios the cold dark matter has a negligible ini- 
tial velocity dispersion (as opposed to "hot" dark matter 
scenarios). This allows us to fully define the particles by 
their initial Lagrangian coordinate q with a unique initial 
peculiar velocity v(q). Then, the equation of motion for 
each particle reads as, once the homogeneous expansion 
of the Universe has been taken into account, 



<9 2 x(q) <9x(q 



ch 



n- 



(2) 



where r = J At /a is the conformal time (and a the scale 
factor) and H. = dlna/dr the conformal expansion rate. 



The gravitational potential <f> is given by Poisson's equa- 
tion 



(3) 



where fi m is the matter density cosmological parameter 
and (5(q) = (p — ~p)/p the matter density contrast. It is 
to be noted that in this expression the Laplacian is taken 
with respect to the x coordinates while the fields are nat- 
urally given as a function of q through the expression of 
the displacement field. We assume here that the den- 
sity contrasts vanish at initial time. The conservation of 
matter then implies that 



1 + %) 



1 



J(q) 



with J(q) 



det 



<9x 
9q 



(4) 



Then, by taking the divergence with respect to the Eule- 
rian coordinate x of the equation of motion ([2]) we obtain 



J(q)V x . 



<9t 2 8t 



= ^ m W 2 (J(q) - 1) 



(5) 

where we used Poisson's equation. As in the Eulerian 
case, it is convenient to introduce the time coordinate r\ 
and the function /(r) defined from the linear growth rate 
D+(t) as, 



r] = In D + (t) 



d In D, 



dlnL> 



+ 



(6) 



d In a TCAt 
The linear growth rate D+(t) is the growing solution of 



A 2 D^ 
dr 2 



H- 



At 



(7) 



which we normalize as -D+o 
reads as 



1 today. Then, Eq.© 



J(q)V, 



2/ 



2 (J(q)-l) 
(8) 

where we note with a prime the partial derivative with 
respect to time rj. 

In the following we will restrict the calculations to the 
Einstein-de Sitter case for which f2 m // 2 = 1. It is to be 
noted however that for all models of cosmological inter- 
est we have f2 m // 2 — 1 so that this assumption is very 
mildly restrictive Q. Thus, up to a good approxima- 
tion, our results can be extended to ACDM cosmologies 
by substituting for the appropriate linear growth rate 
D + (t). Then, the dependence on the cosmological pa- 
rameters is fully contained in the time-redshift relation 
7](z). 

Equation {SI can be written in matrix form as 



Tr 



com(g) 



<9*" 
<9q 



Id*' 
2~dq~ 



= |(J(q) - 1) 



(9) 
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where com(|^) is the comatrix of (dx/dq). It is also 
given by: 



Thus, Eq.© is the form of the equation of motion (J2|) 
written in terms of the Lagrangian displacement field W 
alone. However, it is not sufficient to fully determine 
the dynamics as can be noticed from the fact that we 
only used the potential part of Eq. © when we took the 
divergence in Eq. (J5J) . Thus, we must supplement Eq.© 
with the rotational part: 



d 2 ^(q) 
<9r 2 



H 



0*(q) 



= 0. 



(11) 



As is well-known from the Eulerian perturbation theory, 
the rotational part of the Eulerian peculiar velocity field 
v decays in the linear regime and a curl-free initial veloc- 
ity field remains potential to any order in perturbation 
theory [l|, Q (bu t vorticity will be generated by shell- 
crossings, see [TtJ for an estimation of this effect). Then, 
one usually restricts the dynamics to the case of irrota- 
tional initial velocity fields, V x x v = 0, so that Eq. lfnj) 
simplifies to: 

hence « _ ffiW, (12) 



dxi 



which is of first order over time. In matrix form this 
constraint implies that [l8l |. 



dq 



com(-) 



is a symmetric matrix. 



(13) 

In three-dimensional space Eqs. (fT3|) are cubic in and 
W (in general they are of the order of the number of space 
dimensions). However, it is possible to derive equivalent 
equations that are quadratic in W whatever the dimen- 
sionality of space. They can be obtained through the 
introduction of the velocity potential, T, which the ve- 
locity field is assumed to derive from, ^(q) = dT/dxi, in 
x coordinates. Expressing T in term of ^Sf and imposing 
that d 2 T I dqidqj is symmetric leads to an equivalent set 
of equations of lower order in ip, [l9| |. These equations 
can also be derived explicitly from Eq. (|12p by multiplying 
it by (dxi/dq m )(dxj/dqi), 



dx t <9^(q) _ dxj dtfj-(q) 



dq m dq e dq £ dq rn 
a constraint that in matrix form states that, 

T 

/ (> x \> (CI 1 \ 

is a symmetric matrix. 



Oq 



dq 



(14) 



(15) 



Equations (Ti^ -iriS } are quadratic over \P hence they are 
more convenient to use than Eqs . (H2J)- (|T3J) in three (or 
more) dimensions 21 1. 



B. Linear regime 

The first stages of the dynamics take place at a time 
when the deviations from the Hubble flow are small. 
Then, the equations of motion can be linearized over the 
displacement field 'J. From Eq.(j4]) the Jacobian J(q) 
then reads up to linear order, 



J L (q) = 1 + Tr 



dq 



1 + = 1 



(16) 

where we note with a subscript L all linear quantities. 
Note also that hereafter we define as the partial 
derivative of the displacement field with respect to La- 
grangian coordinates, 



and we introduced its divergence —re, 



«(q) = -V q .*(q)=-X] 



5* 4 (q) 



(17) 



(18) 



It is to be noted that at linear order k is nothing but the 
density contrast. Its time derivative is proportional to 
the velocity divergence. The motion equation © natu- 
rally reads at linear order, 



1 3 
«i(q) + o K i,(q) = o K L(q), 



(19) 



where we recover the two well-known growing and decay- 
ing linear modes: 



k+ = e n and k_ 



-3t7/2 



(20) 



In the following, we shall assume that the initial con- 
ditions are such that only the linear growing mode is 
present (but it would be possible to set different initial 
conditions): 

«z(q>?7) = e n K (q) hence S L (q, 77) = e r '/t (q). (21) 

Note then that at this order the constraint, Eg . (TIB")) . 
implies that ^' L (q) is curl- free in q coordinates, V q x 
^z,(q) = 0; an d so is the linear displacement field. It is 
then entirely determined by its divergence k. 



C. Correlators and propagators 

Because of the mathematical structure of the theory, it 
is obviously very convenient to rewrite the motion equa- 
tions in Fourier space. The Fourier components of the 
field are defined as, 



«(k) 



d"q 

(2tt)" 



ik.q k . 



«(q), 



(22) 
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where d"q is the n— dimensional volume element. The 
Fourier components of the linear displacement field can 
be easily written in terms of the Fourier modes of the 
divergence field, 



* L (k) = i^Mk). 



(23) 



Note that because of the assumed statistical homo- 
geneity and isotropy of space, ensemble average of prod- 
ucts of two Fourier modes vanish for modes that do not 
sum to zero. This property holds for equal as well as un- 
equal time correlators. In the following, we furthermore 
consider Gaussian initial conditions. As it will turn out, 
this is a crucial property. It indeed determines the dia- 
grammatic structure and the contributions to the quan- 
tities of interest. Within this assumption the entire sta- 
tistical properties of the initial density field are defined 
by its power spectrum, Po(k), such that: 



(K (ki)«;o(k 2 )) = (5 D (ki + k 2 )P (fci), 



(24) 



where ( . ) represents ensemble averages over the statisti- 
cal process at the origin of the large-scale structure. 

If the notion of power spectrum has been widely used 
in theoretical and observational cosmology since the early 
eighties, the notion of propagator is relatively new. It has 
been introduced in Q (see also [2(| for the more general 
notion of response functions). By definition it represents 
the ensemble average of the functional derivative of a 
given cosmic field component with respect to the initial 
field value. What we will be interested in here is the 
propagator between an initial convergence mode fto(k) 
and the final convergence mode n{k! 1 rj) (the one defined 
with respect to the rotational parts vanishes for parity 
reasons in case of rotational- free initial conditions). As 
fc(k' ,rj) is the result of a complex nonlinear process, it 
is formally a functional of the whole set of the initial 
density modes (only in the linear regime does it only 
depend on the same k mode). We can then introduce 
the functional derivative of n(k' , 77) with respect to «o(k): 
9rt(k',?7)/9«;o(k). This is a stochastic quantity whose 
ensemble avera ge d oes not vanish for k = k'. It defines 
the propagator [22|, 



= S u (k-k')G(k,r ] ). 



(25) 



The goal of this paper is precisely to investigate the be- 
havior of the propagator G(k,rj). In the linear regime 
the functional k[kq] is trivial and given by Eq. (|21[) which 
implies that, 



Gl(M) 



(26) 



From a perturbation theory point of view, the functional 
k[kq] can be expanded in terms of the initial convergence 
field, 

K(k,r])=jr f d n wi...d n w p <y D (k-^Tw, J 

P =iV V 1=1 J 

x f (?) (wi, . . . , w„;ry) k (wi) . . . k (w p ) (27) 



where the kernels J-^ are symmetric functions of wave 
modes. They are determined by the motion equations for 
k and lo. Then we have 



G(k, V ) = J2 /d n wi...d n w p _i 



xpf^fwi, . . . , Wp-tjk; 77) (k (wi) . . . Ko(wp-i)) (28) 

(the ensemble average of the r.h.s. of this equation en- 
sures that Yli=i w i = so that fc(k— Yli=i w ») ^ s trans- 
formed into Srj(k — k')). Such an expansion can be rep- 
resented in a diagrammatic way by taking advantage of 
the Gaussian initial conditions. This can serve as a basis 
for resummation schemes. We shall illustrate this con- 
struction for the 2D Lagrangian dynamics first. 



III. 2D DYNAMICS 

The aim of this section is to derive explicitly the mo- 
tion equations for the Lagrangian 2D dynamics and to 
explore the resulting propagator properties. Since its 
mathematical structure is simpler than for the 3D case, 
it serves to illustrate the method we develop here to com- 
pute the propagators. 



Decomposition over curl-free and 
divergence-less parts 



We investigate in this section the simpler case of a 
two-dimensional dynamics. This corresponds to pertur- 
bations with ^3 = that do not depend on the third 
coordinate, qs or X3. Therefore, the nonlinear dynamics 
is restricted to the plane (ei, e 2 ) and particles exactly fol- 
low the Hubble expansion along the third axis 83. Then, 
it is convenient to decompose the Lagrangian displace- 
ment field over a curl-free part x an d a divergence-less 
part A as 




9x I d\ 

dqi ~ dqi 
dx d\ 
dq 2 dqi 





- V q .x + V q x (Ae 3 ). 

(29) 

Here and in the following we note x the 3-dimensional 
vector product. Then, the divergence —k reads 

« = -V£x, *(k)=fc 2 X (k). (30) 

In a similar fashion, we define the vorticity as 

oj = -V q A, w(k) = fc 2 A(k). (31) 
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Then, the equation of motion (J9j) reads in Fourier space 
as 

13/" 
«" + -«'--« = / dkldk2 M kl+k2 ~ k ) 

x ^a(ki,k 2 ) 



13 13 

Kl(«2 + o K 2-T K 2)+Wi(w2+ ~W 2 --W 2 ) 



+/3(ki,k 2 ; 



1 13 

ui(k%+-k' 2 )-Ki(u)%+-uj' 2 ) + -kiui2 



(32) 



where we noted m = «;(kj),a;i = w(kj), and we intro- 
duced the symmetric kernels 



a(ki,k 2 ) 



/3(ki,k 2 ) 



det(k 1; k 2 ) 2 
(k 1 .k 2 )dct(k 1 ,k 2 ) 

"1 "'2 



(33) 
(34) 



with 



det(ki,k 2 ) = fci,ifc 2 , 2 - fci, 2 fc 2 ,i = e 3 .(ki x k 2 ). (35) 

It is to be noted that, unlike their Eulcrian counterparts, 
these kernels only depend on the relative angle between 
the wave modes. 

Equation (|32|) can be written in integral form by using 
the Green's function G{rj,rf) that is solution of 



d 2 Id 3\ n . ,, , . ,, 



It reads as 

2 " 



e in-v') _ e -3(»j-V)/2 



(36) 



(37) 



where 6(j) — 7/) is the Heaviside factor which enforces 
causality. This constraint fully determines G{r],r)') 
(whereas Eq. (|3"6"|) alone does not select between advanced 
and retarded propagators or combinations of both). Of 
course, in Eq. (|37[) we recognize the two linear modes of 
Ea. (|2H)l . Thus, we can write the solution of Eq. (|3"2")) as: 

n = K L + f dr/Qfai/) f dkxdk 2 5 D (ki+k 2 -k) 

J —OO J 

x Ja(ki,k 2 ) 



13 13 

«l(«2 + 2 K 2-^ K 2) +W 1 (W2+-W 2 - ^Wa) 



+/3(ki,k 2 ) 



1 13 

Wl(K2 + 2 /t 2)-«l(w 2 +-W 2 ) + -KiW 2 



38) 



where all terms in the brackets are taken at time rf in 
the past. 

On the other hand, the curl-free Eulerian velocity con- 
straint (fT"3l) reads as 



w' = J dk 1 dk 2 5 D (k 1 + k 2 - k) {a(k 1 ,k 2 )[Kiw 2 ~ ^1^2] 

+P(k l7 k 2 )[K lK ' 2 +LU 1 U J ' 2 }}. (39) 



From Sect. Ill Bl we can see that the linear vorticity van- 
ishes, wj, = 0, and Eq. (f3"9")) can be integrated as 

u = f dr/ / , dk 1 dk 2 fc(ki+k2-k){a(k 1 ,k 2 ) 
x [kiuj' 2 — ljik 2 ] + /3(ki, k 2 ) [kik' 2 + uilu 2 }} . (40) 



In Eqs. (|3"8")) and iPTO)) we have set up the initial conditions 
at time r]i —* — 00. It would be possible to keep r/i fi- 
nite, but this introduces extra terms in the perturbative 
series for k and u> that involve the decaying mode k_ of 
Eq.®. By contrast, from Eqs. ([55]). ([30]). the nonlinear 
quantities k and 10 can be written as a perturbative se- 
ries over powers of the linear growing mode e v Kq , such 
that the term of order p factorizes as e pv K^'(k), as in 
the standard perturbation theory. 

The kernels a(ki,k 2 ) and /?(ki,k 2 ) obey the symme- 
tries 

a(k 1 ,k a ) = a(k a ,k 1 ), ^(k 1) k 2 ) = -/3(k 2) k 1 ), (41) 



as seen from Eqs. (l3"3")) -(f3"5' j) . This is consistent with the 
fact that k and \ are scalars whereas to and A are pseu- 
doscalars, as seen from Eq. (|2"9l (so that V q x (Ae 3 ) is a 
vector like VP). Then, under parity V we have: 



V : k 



Av, id — > —Ld, a 



a, (5-* -p. 



(42) 



Equation (|4"2"|) actually implies that the kernels a,/3, sat- 
isfy Eq. (|4Tj) , since the exchange of basis vectors ei «-> e 2 
can be written as a rotation followed by a reflection. 
Then, the symmetry (|42l) directly determines which ker- 
nel a or (3 is associated with a factor such as kk or kw in 
Eqs.dSU) and 



B. Diagrammatic representation 



The equations (f3"51 |4"0")) have a simple diagrammatic 
representation which illustrates the fact that the func- 
tions are obtained from successive quadratic inter- 
actions. A diagrammatic expansion of equations (|3"8"ll4"0|) 
is presented in Fig. [T] Each open circle stands for a 
linear growing mode doublet {kl,lul} — {e Vj Ko(k.j), 0}, 
whereas the vertex points represent the interaction op- 
erators that can be read out from Eqs. (|3"51 |4"0")) . For 
instance, its first component (the one that represents 
{ft(ki,J7i), K(k 2 ,r]2)} -> «(k,r?)) is 

7m (k, 77; ki, 771, k 2 , 772) = J dr]'g(ri,r)')5n(k - ki - k 2 ) 



Then, one must integrate over the coordinates (kj,T)j) of 
the incoming modes at each vertex. 

As noted in [l3[ , each of these diagrams exhibits one 
and only one "principal line" : a line that runs from the 
initial time to the final time without crossing a circle. It 
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{K(k,T]),W(M)}= 



K (w,) 

=o + 




FIG. 1: Diagrammatic expression of the expansion of the convergence- vorticity doublet, {«;(k, 77), uj(k, r/)}, in 2D dynamics. See 
text for details. 



G(k,u)= < + 



=0 + <- 



-O + <= 



=0 



o + 



o + 




=0 



-o 



=0 



FIG. 2: Diagrammatic expression of the expansion of the propagator G(k, rf) in 2D dynamics. All the contribution up to 2 
loops are included. Note that the last two rows correspond to loop configurations that do not all connect to the "principal 
line" (shown here as a straight horizontal line). These are the contribution we assume to be subdominant. See text for details. 



is then possible to sort the loop terms with respect to 
the number of vertices that are attached to this line. It 
is then expected that in the high-fc limit the dominant 
contribution comes from the diagrams whose number of 
such vertices are larger (see [l3| for details). As seen in 
[lH and recalled below, this can be justified in a certain 
regime if it is possible to have a large separation of scales. 
In the following, we will restrict our calculations to this 
subset of contributions. For instance, for the terms up 
to two-loop order, they correspond to the first row of 
Fig. [21 It is important to note that these diagrams are 
such that the incoming waves are always in the linear 
regime. In 13], the authors were able to resum these loop 
contributions (basically by properly counting them). In 
a Lagrangian description things are made more difficult 
because of the complex nature of the vertices and it is not 
always possible to obtain an explicit analytical formula 
for this resummed propagator. 

However, as shown in [16j . the propagator defined by 
this partial series of diagrams can be seen as the exact 
propagator of a simpler dynamics (that only gives rise to 
these diagrams). The latter can be derived by linearizing 
the equations of motion in a certain fashion. Then, we 
can compute the propagator G(fc,7y) by solving exactly 
this second dynamics and next performing the average 
over the initial conditions. This can be made numer- 
ically without performing diagrammatic resummations. 
We first illustrate this alternative method for the 2D- 
Lagrangian dynamics. 

C. High-fc approximation 

1. Resummation of dominant diagrams 

As stated above, the dominant diagrams are expected 
to be those where all incoming lines to the principal path 



are in the linear regime. Following [16| . such a system is 
described by motion equations similar to (|3"5)) and (I40p . 
where in the terms in the right hand sides we replace all 
terms except one by their linear values kl and u>l (the 
latter vanishes here) and sum over all possible choices. 
These equations correspond to a physical system where 
it is legitimate to separate scales, for instance if there 
exists an upper wavenumbcr A so that most of the power 
is associated with small wavenumbers w < A. Then, 
in the limit k 3> A, the evolution of a given mode k 
is governed by the contributions of small wavenumbers, 
k% < A (whence k-i ~ k) or hi < A (whence k% ~ k), in 
the right hand side of Eqs. ([32|) and (|39]) . that are further 
assumed to be in the linear regime. 



The motion equations for the high-k modes then form 
a set of linear equations in presence of a random back- 
ground described by the collection of the low-Wj modes. 
This leaves us with still a complicated system of equa- 
tions to solve. A dramatic simplification can further be 
made because of the high-fc limit. Indeed, since fc 3> Wj, 
Wj denoting the incoming linear wave modes, the wave 
mode k is almost left unchanged along the principal line 
(in other words one is entitled to replace <5D(k' + Wj — k) 
in each vertex point by $D(k' — k)). In this context, 
the motion equations that describe the mode evolution, 
Eqs. (|3li]) and (|3T))) . can by approximated by, 
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K "{Kv) + l; K '(' k ,v) ~t; k (Kv) = J dw k l (w, 77) ja(k, w) 

+/3(k, w 



oj"(k, V ) + -u'(k,r 1 ) 



(44) 



w'CM) = J dw KL (w,7 7 ){a(k,w)K(k,r ? )- W (k,7 ? )]+ / 3(k,w)[ K (k,r;)- K '(k^)]},(45) 

so that high-fc modes now evolve independently on one another. One can easily check that the solution of Eqs. (|4"4"|) - 
(|45p . written as a perturbative series over kq, gives back the principal-path diagrams described above (here with the 
approximation k' = k). 

Of course, we could apply the same procedure to the equations of motion (|38|) .()40 |l . written in the integral form 
[l2[ [l6j. This is equivalent to the differential form used above, but it is less convenient for practical purposes. For 
instance, it is easier to solve numerically the differential equations ([44]) - (|45l) than their integral equivalents which 
require the computation of an integral over all past values to advance to the next time-step. 

Then, we note from Eqs. (|44p - (|45[) that all contributions from the incoming waves Kj,(Wj-) can be factorized out and 
resummed in two distinct bundles of waves, d and (3, defined as, 



d(k) 



such that Eas.(|4i t -B5 |l now read 



1 



dw/tdjwjttf^w), /3(k) = / dw K (w)/3(k, w) 



«"(k, 77) + -n'(k, rj) - -/s(k, V ) = e"d(k) (^"(k, V) + ^'(^ V)) + e"/?(k) (w"(k, 77) + ~u/(k, 77) 

= -eV(k)( K / (k,r ? )- K (k, J 7))+ e "a(k)(a;'(k ) 7 7 )-a;(k, J7 )). 



(46) 

(47) 
(48) 



In other words, the fields «(k) and w(k) depend on 
the linear modes only through the combinations d and 
P. This introduces a dramatic simplification because 
then the ensemble average of Eq. p5|) can be performed 
through a simple average over the two variables a and (3. 
Since Eqs. (j47|) - (|48|) are linear the solution is proportional 
to KqQc). It is convenient to write it as, 



«(k, 77) = e v K (k) k(r)), 
w(k,7y) = e t? K (k) ui(r)). 



(49) 
(50) 



As a consequence, we have 

G{k,ri) = e"G(77) with G(n) = (k(n)). (51) 

We can already note that because a and (3 depend on the 
direction of k only, G(rj) will be completely independent 
of k (since a priori it could only depend on its norm) . 

In the last equation ([ST]) the ensemble average now 
reduces to the computation of the expectation value of 
k(j]) with respect to the distribution of a and f3. We then 
need to explore a bit more the statistical properties of a 
and (3. Using the fact that the linear density field <5z,(q) = 
K^(q) is real, hence kq(w)* = kq(—w), it can be easily 
checked that a and [3 are real numbers. Moreover, we can 
see from Eqs.(j46|) that they are independent Gaussian 
random variables with: 



(a 2 )=3a 2 2 , 2 ) 



'2- 



(&0) - 0, 



(52) 



with, 



dw wPo(w) 



(53) 



Note that 8<r| is also the variance of the density contrast 
(<5(x) 2 ). The joint distribution function of a and (3 is 
then 



V(a,/3) da d/3 



da d(3 
V32ttct| 



■ exp 



6af 



2al 



As a result we simply have, 

G(r))= k(n;a,f3)V(aJ)dadj3, 



(54) 



(55) 



where £(77; a, (3) is the solution of the system (Un- 
written in terms of h and u>, parameterized by the coef- 
ficients d, (3. The calculation of the propagator can take 
advantage of the symmetries (1421) . In particular, we have: 



k(r];a,f3) = £(77; d, -(3), uj{-q;a,f3) = -u>(r);a,-j3). 



We can also note that for /i > 0, 



k(rj; /Lid, /x/3) = k(r] + In d, f3) 
u)(n; /la, fi$) = <2}{r) + ln/x; d, p). 



(56) 

(57) 
(58) 
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2. Behavior of the propagator G(k, rj) 

The asymptotic behavior of G(k,rj) is intimately re- 
lated to the behavior of the solutions of (|47|) - (|48|) for 
finite values of the parameters a and (3. This can be in- 
ferred by inspection of these differential equations. Thus, 
looking for an asymptotic power-law solution, k ~ K^e 1 " 7 
and Cj ~ a) 00 e l " ? , in the limit of large rj where the right 



hand side dominates in Eqs. 
condition 



dH])-([15]), we obtain the 



—(3v av 
= (a 2 +/?>(*, + 1)^+^0, 



which gives the asymptotic modes: 

V\ = 0, h> 2 = — 1, V3 



(59) 



(60) 



In fact, the mode v-i can be removed since Eq. (|47|) can 
be integrated once, as shown in Eq. (|A4|) in the appendix. 
Therefore, when rj 3> 1, k and Cj are expected to be con- 
stant (their value depending on the parameters a and (3 
in a complicated way) because of the mode V\. This im- 
plies that the propagator G obtained from the Gaussian 
integration (|55|) must also be constant at late time. This 
expected behavior assumes that the differential equations 
(|4"Tj) - (|48p do not encounter a singularity at a finite time 
■q. We can check that Eqs. (|47l) - (|48| do not show ex- 
plicit singularities associated with zeros of the coefficient 
of the higher-order terms. Indeed, the determinant of the 
coefficients of highest-order derivatives reads as 



1 



1 



-(3e^ 
-ae r ' 



(1 



ae 



rj\2 



(f3e 



lj\2 



(61) 



which never vanishes if (3 ^ 0. We have checked numer- 
ically that the system of differential equations (|47[) - ([4"5|) 
obeys the behavior described above, with no singularity 
and a constant asymptote a late time. This is depicted 
in Fig. [5] with a 2D plot of k(rj; a, (3) over the plane (a, (3) 
at time 77 = (this 2D plot is sufficient to fully determine 
the behavior of fc(r]; a, (3) thanks to the scaling law (|57p \ 

We show in Fig. [3] our results for the propagator G(rj) 
obtained from the numerical integration of Eq. (|55|) . We 
can see that it first grows until it reaches a maximum 
at rj ~ and next decreases to converge to a constant 
G(rj — 00) ~ 0.8, a behavior qualitatively in agreement 
with the discussion above. Note that at early times the 
rise of G means that the propagator G(k, rj) grows faster 
than the linear prediction (|26[) . until r\ ~ 1. 

The 2D case described above illustrates the power of 
the method based on associating the series of principal- 
path diagrams with a linear dynamics as in Eqs. (|47[) - (|48p . 
Indeed, in the high-Zc limit the dependence on initial con- 
ditions is reduced to a few random parameters (here a 



G(k,j,)/Exp[/7] 




FIG. 3: The propagator G(rj) as a function of rj. It is obtained 
from Eq. (|55[l for a variance 02 = 1. Note that G(rj) depends 
only on the reduced variable rj + In 02 as a consequence of 
15711551) . 



and (3), as can also be read from the diagrams of Figs. Q] 
Then, the ensemble average is reduced to ordinary inte- 
grals (f55|) (instead of path integrals over the field KqQc)) 
and the resummation associated with the infinite series 
of diagrams is obtained by computing the exact solution 
of the differential equations Eqs.([4"7 j) -(|38" ]) , Both steps 
can be performed numerically, as above, since they only 
involve ordinary integrals and differential equations (in- 
stead of functionals of fields). This allows us to compute 
the propagator G(k, rj) even when the diagrammatic se- 
ries cannot be exactly resummed by analytical formu- 
lae (which corresponds to the case when the differential 
equations (I47| - (l48|) have no known explicit solutions). 
Moreover, even in this case, we can obtain exact analyt- 
ical results for the late-time non-perturbative behavior 
of the propagator, as in Eqs. ([60|) . by direct inspection of 
the effective linear equations of motion (|47]) - (j48]) . 

We further discuss the properties of the system ([47]) - 
(|4"5|) in appendix \K\ In particular, we show that taking 
into account the vorticity (i.e. (3^0 and Cj 7^ 0) is nec- 
essary to obtain a well-behaved propagator at late times 
(otherwise a divergence appears), and that the perturba- 
tive series over powers of kq is probably only asymptotic 
(i.e. with zero radius of convergence). 



IV. 3D DYNAMICS 

We now consider the case of the full 3D dynamics. This 
leads to slightly more intricate expressions as we have a 
few more degrees of freedom but we can still follow the 
analysis described in Sect. IIIII for the simpler 2D dynam- 
ics. Moreover, we shall find that the results obtained in 
Sect. Mil remain valid. 

First, as in Eq. (|2"9")) , we can decompose the displace- 
ment field over a curl-free part x an d a divergence-less 
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part A as 



<>\, 



/ dx_ , 9A_ 

' 9gi 9^2 9?3 

d x | gAi 0A 3 

9<?3 9gi 

\ S?3 Oqi dq 2 



Vq-X 



V q X A. 



(62) 



Thus, the rotational part A has now two degrees of free- 
dom: there are three components Ai,A2,A3, but the di- 
vergence of A does not contribute and can be set to zero. 
As in Eqs. (f3"0)) - (f3"l"j) we define the divergence — k and the 
vorticity Co by 



-v qX , 



-V q A. 



(63) 



Then, the spatial derivatives of the displacement field 
read in Fourier space as 



99* 



k-i k j 



= * (J (k) = — ^K(k) - e am ^io m (k), (64) 



kj hi 



dq 



k 2 



where eu m is the Levi-Civita symbol. Then, the equation 
of motion ([9]) reads in Fourier space as 



— K — —K 

2 2 



dkidk 2 fo(ki + k 2 



h 2 h 2 h 2 



1^2 

dkidk 2 dk 3 J D (ki 
1 3 

[(k 2 x k 3 ).(ki x a5i)]« 2 (/«3 + -k 3 - 7K3) + 



fc^-^.k,) 2 „ 1 , 3 



, 2 , 2 [k 2 .(ki x wi)J(/c 2 + -k 2 - ^ K 2) 

1 2 



-%^ Kl [k 1 .(k 2 x(4' + i4))] 

1 2 

det(ki,k 2 ,k 3 ) 



k) 



f det(ki,k 2 ,k 3 ) 2 , „ 1 , 1 



O 1.2 1.2 1.2 



-KlK 2 (,K 3 + -K 3 - ~K 3 ) 



det(ki,k 2 ,k 3 ) „ 1 

, 2 , 2 , 2 «i«2[(ki x k 2 ).(k 3 x (oj 3 + -oj 3 ))\ 

1 



..(65) 



where the dots stand for terms of order ui 2 and uj 3 . We 
do not write these terms here since they will not con- 
tribute to the high-fc approximation. The determinant 
det(ki, k 2 , k 3 ) introduced in Eg. ([63)) is the determinant 
of the 3x3 matrix obtained by putting the coordinates 
of the vectors ki,k 2 and k 3 , in the three columns. It is 
also given by: 



Note that Eq. (|6"5|) is now cubic over hence over k, uj. 
For the constraints associated with the curl- free condition 
(fP2"]) we can use Ea. (fT5)) which is still quadratic. This 
gives: 



det(ki,k 2 ,k 3 ) = (k a x k 2 ).k 3 



(66) 



J 



(k x Co') x k 
fc 2 



dkidk 2 6 D (ki + k 2 - k 



ki x k 2 

p p 



{(ki.k 2 )KiK 2 + Ki[ki.(k 2 x lj' 2 )] + [k 2 .(ki x wi)]4} + ■■ (67) 
I 



where the dots stand for terms of order uj 2 . We can check 
that only the combination kxtJ appears in Eqs.([65 |) -([67 |) . 
Moreover, as in the 2D case where Eqs. (j3"2"|) - (j3"9")l obeyed 
the parity symmetry (|42|) . we can check that Eqs. (f65|) - 
()67|) are consistent with the parity symmetry 



V : k — > k : Co 



■ UJ. 



(68) 



In agreement with Eq. ([62l . A and Co are pseudovectors. 



As in Sect. MI C 1[ the resummation associated with 
principal-path diagrams can be read from Eqs. (|65p - (167|) 
by linearizing over k, Co. This yields 



k" + \ k ' -\ k = J dwe " K o(w) j 

-/ dwduAo (w,« ( U ){ig^£(«» + i«- + |« 



k 2 w 2 — (k.w) 
k 2 w 2 
2 



, 11 1 a (k.w) 



(kx (£" + -£'))]} 



det(k, w, u) r . . ,, 1 _. 3 ., , 
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and 



(k x Co') x k 



= J dwe r, Ko(w) 



w x k 



k 2 J k 2 w 2 

x {(k.w)(V -k) + [w.(k x (uj' - £))]} . (70) 

As for the 2D case, each mode «;(k),u;(k) evolves inde- 
pendently of other high-fc modes and the dependence on 



the initial field ko is reduced to a few random parameters 
that can be written as integrals over n . In order to make 
further progress, it is convenient to write Eqs. (|69| - ((70|) 
in terms of coordinates. Without any loss of generality, 
we can choose k along the axis ej., and Co in the plane 
(e 2 ,e 3 ). Then, Eqs.pt)- ((SI read as 



«" + ~«' - |« = e"(r 22 + 733)(«" + - eV 13 K + \j 2 ) + eV 12 K + \j 3 ) - e 2 "(r 22 r 3 3 - r 2 3 )( K " + \k' + |«) 

13 13 

-e 2?7 (ri 2 r 23 - t 22 t 13 )(lo' 2 ' + ~v' 2 + -u 2 ) + e 2ri {r 13 T 23 - t 33 t 12 )(lo 3 + -u>' 3 + -u) 3 ), (71) 



and: 

J 2 = e ,, T 13 (K'~K) + e v T 33 (uj' 2 ~uj 2 )-e ,, T 23 (uJ 3 ~uj 3 ), (72) 

<4 = -e' ? Ti 2 (K / - re) - e' ) T 23 (w 2 - lo 2 ) + e n T 22 (u 3 - u) 3 ). 

(73) 

Here we introduced the symmetric parameters defined 
by: 



J dw kq( w ) 



U>i Wj 



(74) 



Using the property ko(w)* = Kq(— w), we can see that 
the coefficients are real random numbers. We recover 
the two-dimensional case ([S7 | -(|35 ]l for 

r i3 = 0, uj 2 = 0, a = r 22 , $ = Ti2, (75) 



which behave respectively like spin 0, 1 and 2 com- 
plex numbers with respect to coordinate rotations in the 
(e 2 ,e 3 ) plane. The ensemble average of those quantities 
can be expressed in terms of a 2 defined as, 



with 



8tt f , 2 . , 
— dww P (w), 



(r 2 ) = (Ivl 2 ) = 



'3> 



(82) 



(83) 



t j2 = 0, w 3 = 0, <5 = t 33 , /3 = -n 3 . 



(76) 



Note that there are several symmetry properties. Two 
symmetries extend the property (f5"6"l) obtained for the 
2D case. They read as 



and 



Tl3 



Tl2 



-Ti3j 723 



"Tl2, T 23 



-T23, W 2 



-T23, U> 3 



-^3, 



(77) 



(78) 



where we only write the quantities that change under 
these two symmetries. A further symmetry comes from 
the invariance over a coordinate rotation in the (e 2 ,e3) 
plane. To express it, we can define the following quanti- 
ties, 



T = 



T 2 2 + T 33 



V 

7 = 7e 



V e = T12 + 1T13, 
2i fl _ T22-T33 



1T23, 



(79) 
(80) 
(81) 



while cross-correlations between these quantities vanish. 
We also define the complex vorticity ui as 



UI = -LU 3 + 1UJ 2 , 



(84) 



which is of spin 1 like v. Here we use the fact that, as 
in the 2D case, Eqs. ([7T|l -([73 |l are linear so that we can 
factorize a factor reo(k), as in Eqs. (|4"9")) - (f50")) . Then, the 
reduced quantities k/kq, uJi/kq, are real (since the coef- 
ficients Ty of Eq. (|74|) are real) so that the complex vor- 
ticity (|84|) fully determines the doublet {oj 2 ,uj 3 }. Then, 
the two equations f?2" ]) -(|73| ) can be gathered into 



/ = p0 



e"v(re' - re) + eJ>T{J - ui) + e"7(tJ , - uj)* , (85) 



whereas Eq. ([7T|) reads as 
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-k' - -k = e'"2r( K " + -k') - e" — (u" + -u>') - e^-(uj" + -lj'Y - e 2 "(r 2 - ffMn" + -k' + -k) 
2 2 v 2 2 v 2 ^ 2 ^ 2 ; v A 2 2 y 



i 



-TV, „ 1 . 3 . „ 



(86) 



G(k,^)/Exp[/,] 




FIG. 4: The propagator G?(r;) as a function of r/ for the 3D 
dynamics. It is obtained for a variance as = 1. 



We can see that all terms in Eq. (j85|) are of spin 1, whereas 
all terms in Eq. (f8"6"]) are of spin 0. This clearly shows 
that these equations are invariant through rotations in 
the (e2,es) plane. Moreover, we can check that both 
sides in Eq. (|86|) are real. Obviously, the results depend 
only on the angle difference 8 V — 6* 7 . 

Finally, the scaling laws ([57| -(|58 p also extend to the 
3D case as 



fi > 



rj — In fj,. 



(87) 



As for the 2D case analized in Scct. lIII C 21 we can look 
for singularities associated with zeros of the determinant 
of the coefficients of higher-order derivatives. This gives 
from Eas.([8"51) 



A = [(1 - e"r) 2 - e 2 "|7| 2 r + |e"v + e 2 "(v*7 - rv)| 

(88) 

The determinant A can only vanish if \ j\ = 0, or |v| = 0, 
or 9 V — 8 -y — nir/2 with n integer, which is a region of 
zero measure in the space spanned by the coefficients t, v, 
and 7. We have checked numerically that the differential 
system is otherwise well-behaved and the ensemble aver- 
ages lead to well-defined quantities. As in Eqs. ((59|) - ([6ll)) . 
the asymptotic behavior of the solutions k, 1^2,^3, can be 
read from the differential equations (f7"l"j) - ([73"|) by looking 
for asymptotic power-laws. This yields for the reduced 
variables k, &2,&3 , defined as in Eqs. p5|) - (|5U1) . the three 
asymptotic modes: 



v\ = 0, v 2 



iV23 



-, ^3 



-5 + iV23 



(89) 



Therefore, the reduced propagator G(rj) must go to a 
constant at late times, as for the 2D case. Our numerical 



results are shown in Fig. 2] and we can see that they 
agree with this analysis. Thus, it appears that the 3D 
propagator exhibits the same features as the 2D case, 
with an early rise that is faster than the linear prediction 
and a late-time behavior that follows the linear power- 
law G(rf) ~ e 1 ' (with a "renormalized" amplitude that is 
smaller than unity). 

As for the 2D case (see Eq. (|53"]) ). the key quantity <r| 
that measures the amplitude of the fluctuations and the 
state of gravitational clustering, see Eqs. ([52]) - ([8^)l . is pro- 
portional to the variance of the density field (i5(x) 2 ). We 
can note that for a CDM power spectrum it shows a log- 
arithmic UV divergence (since Pq(w) ~ w~ 3 at high w). 
Therefore, our results rigorously apply to linear power 
spectra with a high-fc cutoff such that erf is finite. How- 
ever, since the fluid description does not hold beyond 
shell-crossing it could be argued that integrals such as 
(|82|) should be cut at the scale associated with the tran- 
sition to nonlinearity in any case. On the other hand, 
within the high-fc approximation studied in this article, 
the quantity 03 of Eq. ([82| should be interpreted as the 
variance of the larger-scale density contrast, rather than 
the variance of the one-point density contrast. Indeed, 
we can see from Eq. ([71)) that the quantity which gov- 
erns the coefficients is the density contrast at the ori- 
gin 5(x = 0) (discarding the angular dependence associ- 
ated with WiWj/w 2 ). This in turn gives rise to Eq. (|82|) . 
Mathematically, the specific role played by the origin is 
related to the breakdown of the invariance through trans- 
lations entailed by the approximation <Srj(ki + k2 — k) ~ 
<5o(ki — k) discussed in section [ill C II However, it is 
clear that within this approximation, based on a sepa- 
ration of scales between low wavenumbers w < A and 
high wavenumbers k 3> A, any point located at a dis- 
tance below 1/ A from the origin could as well be chosen 
as a reference. In other words, within this high-fc approx- 
imation, erf should be understood as the variance of the 
larger-scale density contrast, associated with wavenum- 
bers w < A (and A < fc). Then, in Eq. (f8"2")) we relaxed 
the cutoff A, which is valid for linear power spectra with 
small high-fc power so that the integral converges (and 
the high-fc approximation discussed in section fill C II can 
make sense) . We can see that CDM power spectra are at 
the limit of applicability of this approximation. 

We can note that the same features apply to the Eu- 
lerian description, except that instead of the larger-scale 
density contrast the key quantity is the larger-scale veloc- 
ity. Then, it happens that CDM power spectra are fully 
within the range where the velocity integral analogous to 
Eq. (|8"2"l) converges. 
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V. DISCUSSIONS 

We have applied to the Lagrangian formalism a re- 
summation scheme developed in [13j within the Eulerian 
description. This is based on the resummation of a cer- 
tain type of diagrams, called "principal-path diagrams" 
m [13jj , that may be expected to dominate the dynam- 
ics in a high-fc limit. In the Eulerian case, these dia- 
grams can be explicitly computed, order by order, and 
resummed, as one can recognize the exponential function 
in the series expansion obtained in this manner. This 
leads to a Gaussian decay of the form e _e Vk cr "/ 2 at high 
fc. 

It is more difficult to apply the same method to the 
Lagrangian formalism, as the diagrams have a slightly 
more intricate expression and one cannot identify from 
the series a well-known mathematical function. However, 
as shown in [lg], it is possible to identify this resumma- 
tion with the solution of an effective linear dynamics. 
Then, instead of computing explicitly all diagrams and 
next resumming their contributions (which amounts to 
solve for this effective equation of motion as a pertur- 
bative series), one can directly solve for this simpler dy- 
namics. In this article, we have applied this technique 
to the Lagrangian description. We have shown that it 
is very powerful as it can be used even when no explicit 
analytical solutions can be found (but one can still solve 
numerically the relevant differential equations). More- 
over, even in such cases, it is possible to obtain the exact 
exponents (as defined by this partial resummation) of the 
late-time regime, by looking for the asymptotic modes of 
the linear differential equations. Then, we have found a 
late-time power-law behavior for the propagator, which 
actually simply follows the linear growth e v albeit with 
a "renormalized" amplitude slightly smaller than unity. 
This is quite different from the Gaussian decay obtained 
in the Eulerian case. 

For comparison, let us briefly recall how this method 
applies to the Eulerian case [l6|. In this case, the solution 
to the effective linear dynamics can be derived explicitly 
and it reads as 

d(k lV ) = e v S (k)e eV&E(u \ (90) 

with 

a E (k) = [ d"w^<5 (w), (91) 
J w z 

(using notations that straightforwardly extend those used 
throughout the paper). For Gaussian initial conditions 
the ensemble average of this expression can be easily com- 
puted. It leads to the following propagator, 

G E {k,ri) = e "e e2 " <<iE(k)2)/2 . (92) 

Because of the k-dependence of a.g;(k), one obtains a 
Gaussian damping of the form e~ e <T « / at high fc, 
with cr 2 = 1/n J (d n w / w 2 )Pq(w) (n is here the number of 



space dimensions). As discussed above and shown in de- 
tails in previous sections, our calculations in Lagrangian 
space do not give a closed form for the propagator but 
allow nonetheless to describe its properties exhaustively. 

Eulerian and Lagrangian calculations prove to lead to 
quantitatively very different results. Whereas the decay 
found for the Eulerian case exhibits a Gaussian tail with 
a strong fc-dependence, in Lagrangian variables the prop- 
agators are essentially fc-independent with no significant 
decay at late time. After a stage of accelerated growth, 
followed by a transitory slow-down, the high-fc modes 
growth is indeed found to be simply slightly retarded 
and still growing as e'' as the linear growth rate. The sit- 
uation is the same in 2D and 3D cases. The delay is only 
slightly less important for the 3D case. The independence 
on wavenumber fc in the Lagrangian case directly follows 
from the fact that the kernels a and (3 of Eqs. (j3"3"|) - ([34]) . 
that appear in the 2D equations of motion (|3"2"|) and (|39p , 
are homogeneous functions of their two arguments ki and 
\i2- they only depend on relative angles. This also holds 
for the 3D dynamics, as can be checked in Eas. (p5|) - ([S7|) . 
Therefore, this property is not restricted to the partial 
resummation associated with "principal-path" diagrams. 
In a similar fashion, the dependence on fc obtained in the 
Eulerian case is due to the non-homogeneous character 
of the kernels a and /3 that appear in this framework, 
which can also be seen in Eq. (f9"Tj) . 

As seen in the previous sections, another distinctive 
feature of the Lagrangian description is the important 
role played by parity symmetries. Indeed, whereas in 
the Eulerian framework the two quantities of interest, 
the density and the velocity divergence, are true scalars, 
in the Lagrangian framework we must take into account 
both curl-free and rotational parts of the displacement 
field (in Lagrangian space q), as a curl- free Eulerian ve- 
locity field does not translate into a Lagrangian curl-frcc 
displacement field beyond second order. We have shown 
that keeping track of the vorticity degrees of freedom is 
necessary to obtain a well-defined propagator in the non- 
linear regime. 

How to reconcile these results? Although Eulerian and 
Lagrangian descriptions are ultimately equivalent, the 
objects we have computed are clearly distinct. In the 
nonlinear regime modes in Lagrangian space cannot be 
directly mapped to those in Eulerian space. One should 
then not be too surprised to find quantitatively different 
results. What we have computed here is in essence the 
leading effect of a random background on the growth of 
structures, assuming scales can be well separated (e.g. 
that the wavelength of the background modes are much 
larger than the modes of interest). It turns out that 
the Eulerian modes are sensitive to the large-scale dis- 
placement field at leading order, whereas the Lagrangian 
modes are not. These large-scale displacement fields are 
responsible for the decay of the Eulerian correlators at 
large time separations. Indeed, the modes behave as if 
they were randomly advected by the large-scale displace- 
ments (13 , [ill ] . Basically, everything happens as if small- 
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scale structures were moved around; and because they 
occupy a different location in real space, their correlation 
with the initial field decay. In Lagrangian space, modes 
are not affected by such displacements (by construction, 
the convergence k and the vorticity u> are not sensitive 
to a uniform translation, being related to derivatives of 
the displacement field taken as a function of the initial 
conditions) . They are more directly sensitive to the den- 
sity field. Thus, as discussed in section [iVl whereas the 
Eulcrian propagator is governed by the amplitude of the 
larger-scale velocity, the Lagrangian propagator is gov- 
erned by the amplitude of the larger-scale density. Then, 
the leading effect resembles more a tidal effect. What 
we have found is that modes are not disrupted by the 
accumulation of those tidal effects, at this order of the 
calculation, e.g., the results displayed in Figs. [3] and Q] 
suggest that there is no true loss of memory nor efficient 
relaxation associated with the gravitational dynamics. It 
is not clear then how this loss of memory - which is ex- 
pected to happen eventually in the nonlinear regime - 
could take place. Whether it can be described with the 
help of additional diagrams [23j, from terms beyond the 
high-fc limit; or whether we have to go beyond shell- 
crossing (which breaks the analyticity of the Jacobian) 
to capture possible relaxation effects, is yet unclear. 
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Then, eliminating & from Eqs. (|A2p - (|A4p gives the 
second-order equation for k: 



2(1 -3aD)D 



i d 2 k , 

2 — + { 7-15aD)D- 



aD(l - aD) 
(1 - aD) 2 + 0Df 



[12k - 15] = (A5) 



We can note that for (3 — the divergent part k de- 
couples from the vorticity cu and the solution of Eq. (jAlj) 
can be written as 



k(D; a, 0) = 2^(1, 3/2; 7/2; aD). 



(A6) 



It exhibits a singularity at the point D = 1/a (for a > 0), 
in agreement with Eq. (|6ip . but it actually remains finite 
at this point and has a well-behaved analytic continua- 
tion beyond, as seen in Fig. 

Here we can note that writing the high-fc resumma- 
tion in terms of the differential equations (T4"4"l) - (|4"5)) . and 
the propagator with the integral representation (|55p . is 
a key ingredient to obtain the asymptotic behaviors. In- 
deed, computing G(k, rj) from its diagrammatic expan- 
sion, which amounts to expand the integrand in Eq. (155p 
over powers of a and /3, leads to an asymptotic series 
with zero radius of convergence. For instance, for (3 — 
we directly obtain from Eq. (IA6l) 



k(D;a,0)=J2 
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(2p 



3)(2p + 5) 



(A7) 



which gives after we average over a: 



APPENDIX A: ALTERNATIVE EQUATIONS 
FOR THE 2D DYNAMICS 



(^&m & = ± ( ^uZ. ,^ D^. ( A8) 



^(4p + 3)(4p + 5) 



Here we explore in a more details the properties of 
system (g7])-(|25]). With the change of variable D — e v it 
yields, 



d 2 k 
dD 2 



7 dk 

-D- — = aD 
2 dD 



D< 



d 2 k 



dD 2 



7 dk 
2°dD 



+(3D 



D 



d 2 LJ 

dD 2 



7 d* 
2 dD 



D- 



dui 
'dD 



dk 



duj 



[3D —— + aD —— 



with the initial conditions: 



D -> 



k = 1 



dD 

3 



dD 



aD, 



Equation (|A1|) can be integrated once to give 



D- 



( dk 
'dD 
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dk 


3 




- = aD 
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D dD- 


h 2 K 



D- 



dui 
'dD 



(Al) 
(A2) 

(A3) 



3 
2 

(A4) 



2 W 



This asymptotic series describes the early rise of G but 
it cannot give (without ambiguities) the late-time relax- 
ation to a constant. 

This behavior emerges because of the existence of the 
second degree of freedom associated with the vorticity, 
[3. To take it into account, one may look for a solution of 
Eqs. (IA2p - (|A4p as a perturbative series over powers of D, 
as in Eq. (|A7l) . Then, computing the first few terms or 
looking at simplified cases suggests that the nonzero vari- 
ance of (3 decreases somewhat the coefficients of Eq. (|A8|) 
but they remain positive and fastly growing (it typi- 
cally modifies Eq.(jA8| by changing the factor (3<r|) 2p 
into (2cr|) 2p , because of Eq . ([521 ) . Thus, as expected the 
vorticity slows down the rise of the propagator G(k, rj) 
but its magnitude is not sufficient to make it decay with 
respect to the linear propagator at early times. Never- 
theless, it is necessary to take into account the vorticity 
to obtain the late-time behavior of G(k,rj). 
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FIG. 5: The divergence k(r] = 0; a, 13) (left panel) and the vorticity (right panel) as a function of a, 0. The divergence reaches 
a constant at large radius, \J at 2 + /3 2 — > oo, for a fixed polar angle. The divergence is found to be continuous and even with 
respect to f3; the vorticity is found to be discontinuous along the critical half-line e v a > 1, = 0, and odd with respect to f3. 
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